Electron-vibration coupling constants in positively 

charged fullerene. 

Nicola Manini 1,2 * Andrea Dal Corso 1 ' 2 } 
Michele Fabrizio 1 ' 2 ' 3 ! and Erio Tosatti 1 ' 2 ' 3 ^ 
Istituto Nazionale di Fisica della Materia (INFM), 
Unitd Trieste SISSA, 
2 International School for Advanced Studies (SISSA), 

Via Beirut 4, 1-34013 Trieste, Italy 
3 International Centre for Theoretical Physics (ICTP), 
P.O. Box 586, 1-34014 Trieste, Italy 

April 10, 2001 



Abstract 

Recent experiments have shown that Cgo can be positively field-doped. In 
that state, fullerene exhibits a higher resistivity and a higher superconduct- 
ing temperature than the corresponding negatively doped state. A strong in- 
tramolecular hole-phonon coupling, connected with the Jahn- Teller effect of 
the isolated positive ion, is expected to be important for both properties, but 
the actual coupling strengths are so far unknown. Based on density functional 
calculations, we determine the linear couplings of the two a g , six g g , and eight 
h g vibrational modes to the H u HOMO level of the Cqq molecule. The cou- 
plings predict a D5 distortion, and an H u vibronic ground state for C^ . They 
are also used to generate the dimensionless coupling constant A which controls 
the superconductivity and the phonon contribution to the electrical resistivity 
in the crystalline phase. We find that A is 1.4 times larger in positively-charged 
Ceo than in the negatively-doped case. These results are discussed in the con- 
text of the available transport data and superconducting temperatures. The 
role of higher orbital degeneracy in superconductivity is also addressed. 
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1 Introduction 



Recent breakthrough experiments (Schon et al. 2000) have shown that a monolayer 
of C 60 can be positively field-doped, a goal hardly realized chemically so far (Datars 
and Ummat 1995). In that state, fullerene exhibits a higher resistivity than for 
negative doping, and becomes superconducting with critical temperatures that can 
exceed 50 K, about a factor 5 higher than the corresponding negative field-doped 
state. The general belief is that superconductivity in the fullerenes is related to 
a strong intra-molecular electron-phonon coupling, connected with the Jahn- Teller 
effect of the isolated ion (Antropov et al. 1993). Unlike the negative doping case, 
where both calculations (Antropov et al. 1993, Varma et al. 1991, Lannoo et al. 
1991) and fits to data (Gunnarsson et al. 1995) exist, no quantitative evaluation of 
the actual Jahn- Teller coupling strengths is so far available for the positive fullerene 
ions. 

In this work we undertake the task of determining the electron-vibration linear 
couplings for the H u HOMO level of the Cgo molecule, along with that of spelling 
out some of the consequences for resistivity and for superconductivity. For that 
purpose we use density functional electronic structure calculations, yielding accurate 
molecular vibration frequencies and eigenvectors for the two a g , six g g , and eight h g 
modes that couple linearly to the electronic H u state. Knowing the form of the Jahn- 
Teller coupling matrices, we distort the molecule and extract the coupling constants 
from the calculated level shifts and splittings. As a parallel check, we repeat a similar 
calculation for the negative, electron doping case, where the a g and h g modes couple 
linearly to the T lu LUMO level of the Cqq molecule. The couplings obtained for 
negatively doped C 60 are rather similar to those that can be found in the literature, 
and just represent a fresher, state-of-the-art theoretical determination. The hole- 
vibration couplings of positively doped Ceo are new, and can be put to use in a 
variety of manners, including predicting or explaining properties of molecular ions, 
such as photoemission (Briihwiler et al. 1997) and IR/Raman spectra. That is a task 
that we propose to consider in the near future. 

The couplings obtained can also be used to determine the dimensionless electron- 
phonon coupling constant A relevant for the superconductivity as well as for the 
vibron contribution to the high temperature T-linear resistivity in the crystalline 
phase. Comparing values for positive and negative doping we find that for positive 
doping A is a factor 1.4 larger than for negative doping. These results provide a 
starting point for a discussion and comparison with the experimental findings. 

This paper is organized as follows: the notation is set up in Sect. ||; the ab-initio 
calculation and results for the molecular ion are described is Sect. ^ in Sect. | we 
sketch the calculation of the resistivity in a band-degenerate case; Sect. |^ presents 
a formulation for superconductivity in that case; Sect. |B| contains comparisons and 
discussion of experimental data. 
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2 The Jahn- Teller Hamiltonian 



Several theoretical papers (Ceulemans and Fowler 1990, De Los Rios et al. 1996, 
Moateei al. 1996, Moate et al. 1997, Manini and De Los Rios 2000) formulate the 
H <g) (a + g + h) Jahn- Teller (JT) and dynamical JT problem - describing a hole 
in molecular Cqq - though with different notations and conventions. For ease of 
comparison, it is therefore useful to set up explicitly the conventions we use in the 
present calculation. 

The analogies of the icosahedral H © (a + g + h) JT coupling to a spherical T>^ © 
(x>(°+) + 2}( 2+ ) + D( 4 +)^j model were exploited in earlier work, where among other 
things it was shown that surprisingly, and depending on numbers, this dynamical JT 
problem may or may not possess a Berry phase (De Los Rios et al. 1996, De Los 
Rios and Manini 1997, Manini and De Los Rios 1998). The representations T>( L± ^ of 
0(3) map into representations of the icosahedral group X^ as follows (Altmann and 
Herzig 1994): — > a g/u , — > t lg/u , — > h g/u , — > t 2g/u ®h g/u , 

2}( 4± ) — ► g g / u ®h g / U) and so on. This means that the decomposition of the symmetric 
part of the tensor product 

{V^ © £> (2_) } s = V (0+ ^ © V (2+ ^ © (1) 

becomes in icosahedral language 

{K ® K} s = a g ®h g ® (g g © h g ) . (2) 

a g and g g appear only once, while the h g representation modes appears twice in this 
decomposition: one represents a genuine quadrupolar T>^ 2+S) state, while the other 
one derives from a £>( 4+ ) representation. Though this 0(3) picture is suggestive 
(Ceulemans et al. 1994), clearly a quantitative description of Cgg 1 " ions had better 
involve icosahedral symmetry from the beginning. 

In the icosahedral group indeed in the h u ®h u tensor product, the h g representation 
appears twice. This reflects the non-simple reducibility of the icosahedral symmetry 
group. Accordingly, Butler 1981 provides two independent sets of Clebsch-Gordan 
(CG) coefficients 

fc C£M = (H,n;H,u\h,m)M (3) 

which couple an H electronic state (quadratically) with an h vibrational mode (lin- 
early) to give a scalar. Each set of coefficients is identified by a multiplicity index 
r = 1, 2. Since the two h states labeled r = 1, 2 are symmetry- wise indistinguishable, 
the choice of these two sets of coefficients is perfectly arbitrary, as long as they are 
kept orthogonal to each other. This arbitrariness is the source of the different nota- 
tions taken in the literature of this field. Here, we stick to Butler's choice (Butler 
1981), which is basically equivalent to Ceulemans' convention (Fowler and Ceulemans 
1985). Also, we label the states within a degenerate multiplet by the labels of the 
subgroup chain Xh D D 5 D O5. For brevity, we indicate only the O5 index m (m = 
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for a g , m = ±1, ±2 for g g and to = —2, . . . , 2 for h g states) in the labeling of states 
since, for the representations relevant to our problem, the D 5 label is just the absolute 
value of to. 

Given the tabulated CG coefficients, it is necessary for generality to consider a 
linear combination 

h C™ u (a) = cos a h C^ + sin a h C™® (4) 

of the two sets. The coefficient h C™ u [a] coincides with Butler's r = 1 and r = 2 
values for a = and a — | respectively. Different values of a can be compared with 

the conventions of previous authors. For example, a = — arctan ~ —53.3° is 

the case studied by De Los Rios et al. 1996 (where the X n CG coefficient becomes 
equivalent to the spherical (2, /i; 2, u\2, to)); a = ir/2 by Moate et al. 1996; and 
a = 0, 7r/2 by Ceulemans and Fowler 1990, where these cases are indicated as h b and 
h a respectively. The a-dependence of these CG coefficients indicates that - unlike, for 
example, cubic symmetry - belonging to the h g group representation in icosahedral 
symmetry does not determine completely the form of the JT coupling. The mixing 
angle a is also needed for that. In the present case of fullerene, each h g distortion 
mode is thus characterized not only by its frequency and scalar coupling, but also by 
its specific mixing angle — n/2 < a < it/ 2. 

The basic linear Jahn- Teller Hamiltonian for the H u ® (a g + g g + h g ) model is 
conveniently divided into: 

^modcs ( T ) 

Q Ti ) + H^_ v (g ri hu Ti Qn)\ ■ (5) 

T=a, g ,g a ,h g i 

The first term describes the linearly-coupled vibrations in the harmonic approxima- 
tion, 

HLJ^, P,Q) = ^f E(^ + Q 2 J . (6) 

Z rn 

while the second term is the linear coupling to each mode: 

Hl_ v (ghu;, a,Q) = 9 -^Y,Qm c\c. v T C™ U (a) . (7) 

Here, of course, the a dependence is relevant only for the r = h g case. The dis- 
tortion coordinates Q T i m (with conjugate momentum P r j m ) are dimensionless, being 
measured in units of Xo(co T i) = (%/mc uj rl ) l l 2 . The operator creates an electron in 
orbital //(= -2, . . . , 2) of the HOMO H u shell. 

Naturally this form of the coupling Hamiltonian is such that each term represents 
pertinent irreducible representation combinations that are totally symmetrical, i.e., 
scalars, under the icosahedral group. For future applications, it will provide a conve- 
nient form both for perturbative calculations (small g values) and as a starting point 
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for numerical diagonalization methods, such as the Lanczos technique. Here, how- 
ever, we restrict ourselves to a study of the classical molecular distortions. For that, 
it is more convenient (Manini and De Los Rios 2000) to switch to a real representation 
of the vibrational degrees of freedom and orbitals. 

To that end, we apply two standard unitary transformations (Manini and De Los 
Rios 2000), one to the electronic and the other to the vibrational degrees of freedom. 
We define a new set of electronic operators, d m (and consequently their Hermitian 
conjugates d^), as 



Co 



d 
1 
72 



(8) 



dr, 



m 



1,2 



This transformation leaves unchanged the (diagonal) coupling to the a g modes, which 
takes the final form: 



(9) 



fj, V 



with 



bav and q = Qo- 



The second (similar) transformation is applied to the vibrational coordinates of the 
g g and h g modes: 

(10) 



Qo = qo (h g modes only) 

Qm \ = (~l) m (l i \ ( q m 



m 



1,2). 



The harmonic part is left unchanged by this transformation, while the interaction is 
transformed into 



H^ v (ghuj, a,q} = ^ Y, Qm d l d ^ V ^ (m) («) ( r = 9g, h g) ■ 



(11) 



The 5x5 coupling matrices V r ^ are combinations of the CG coefficients. Their 
explicit expressions [we use the shorthand s for v^3 and omit the explicit indication 
of dependence \ hg ( m \a)\ are the following: 
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In the static JT effect, the kinetic term in P? im is ignored, and the problem is to 
study the five Born-Oppenheimer (BO) potential sheets given by the sum of each 
eigenvalue of the electronic problem plus the harmonic restoring forces. The H u ® 
a g part has a purely diagonal coupling matrix. As it does not split the electronic 
degeneracy, it is trivially separated from the coupling to the other modes and can be 
treated as a displaced oscillator. 
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When a single electron (hole) is placed in the H u orbital, the molecule distorts in 
such a way that the lowest (highest) BO sheet is lowered (raised) in energy as much as 
possible. The coupling to a g g mode leads to 10 equivalent absolute minima of D 3 local 
symmetry (Ceulemans and Fowler 1990, Manini and De Los Rios 2000) of the BO 
potential. The optimal distortion is realized for |qf 9 | = g 9g /3, with a corresponding 
potential energy lowering of hu gg /18. The r = 2 part of the coupling to the h g 
modes (corresponding to the sin a terms in V^ 9 ( m ) above) contributes to these same 
minima in an equivalent way, with g g replaced by sin a. However, the r = 1 
component (cos a terms in V^ 9 ( m )) of the coupling favours the six classical stable 
minima of local D$ symmetry (Ceulemans and Fowler 1990, Manini and De Los Rios 
2000). The optimal distortion at these minima is |g^ 9 | = gn g cosa/v^o", for an energy 
lowering of g\ cos 2 a ?ku/10. 

The simultaneous linear coupling to several modes will generally lead to a cumu- 
lative distortion and to an energy gain which is the sum of the individual energy 
gains. However, the form of the coupling (|TT| ) prevents the molecule to gain energy 
through both kinds of couplings. The system shall choose between a D3 and a D$ 
distortion, depending which one is energetically more convenient for given specific 
values of the couplings, vibration frequencies, and H u orbital electronic filling. The 
calculation of the following section determines in particular which one of the two 
types of distortions prevails in C^ . 



3 Calculation of the couplings and results 

We compute the molecular electronic structure within the density functional theory 
(DFT) in the local density approximation. The Ceo molecule is repeated periodically 
in a large fee supercell lattice. The conventional supercell side is a = 18.5 A, so 
that the distance between the centers of two neighboring copies of the molecule is 
13.1 A, suitably much larger than the fully relaxed equilibrium (opposite C-C) ball 
diameter, about 7.053 A. Since we aim at describing the single molecule - and indeed 
our molecules are well isolated - no sampling of the Brillouin zone is called for, and 
calculations of the charge density are done using the k=0 wavefunctions. We use 
ultrasoft pseudopotentials (Vanderbilt 1990) for C (Favot and Dal Corso 1999). The 
plane- waves basis set is cut off at E cnt = 27 Ry (charge density cutoff = 162 Ry). 
Test calculations with higher cutoff or larger cell size a gave equivalent results. 

Based on this electronic structure calculation, we used next density functional 
perturbation theory (Baroni et al. 1987) to compute three independent rows of the 
dynamical matrix. Icosahedral symmetry is then used to recover the full matrix, 
which determines the normal modes £j )S and frequencies u>i (Giannozzi and Baroni 
1994) of the molecule. We obtained frequencies (see Table |l|) in good agreement 
with experiment (Prassides et al. 1991, Zhou et al. 1992), as well as with previous 
calculations (Giannozzi and Baroni 1994, Negri et al. 1988, Kohanoff et al. 1992). 



7 



mode 


Experim. 


Giannozzi & Baroni 


this work 


a g\ l ) 


4yo 


4yo 


rnn 

OUU 


a g {2) 




i F;n/i 
lbU4 


1 CI 1 

Ibl 1 


h (-\ \ 


971 


zoy 


OR1 


h g {2) 


4o I 


/I o c 
4ZD 


/ion 
4zy 


hg{6) 


71 n 
< 1U 


71 1 
(11 


71 S 

/ lo 




77/I 

( (4 


7SQ 

loo 


7fi A 

1 o4 


h n (5) 


1099 


1120 


1119 


h g (6) 


1250 


1281 


1275 


h 9 (7) 


1428 


1450 


1456 


h g (8) 


1575 


1578 


1588 


*i«(l) 


527 


527 


533 


*i«(2) 


576 


586 


588 


fi«(3) 


1183 


1218 


1212 


*i«(4) 


1428 


1462 


1469 



Table 1: Eigenfrequencies (in cm -1 ) for the Raman- and IR-active vibrational modes 
of C 60 molecule: comparisons with experimental values (Prassides et al. 1991, Zhou 
et al. 1992, Gunnarsson et al. 1995) and previous calculation (Giannozzi and Baroni 
1994). 

To determine the e-v couplings for the linearly coupled modes, we proceed subse- 
quently to displace the atomic positions from the equilibrium position along each of 
the normal modes, choosing a suitable eigenvector combination in the linear space 
of each degenerate vibration. In particular for each h g mode we selected the go dis- 
placement, corresponding to the totally symmetric combination £j j0 of the distortions 
£i tS with respect to an (arbitrarily chosen) D$ subgroup of the molecular symmetry 
group. The five initially degenerate (really, only nearly degenerate, owing to a weak 
cubic splitting due to the artificial supercell lattice) H u Kohn-Sham eigenvalues split 
under this distortion with a pattern given by the eigenvalues of ^V h<! (°\ We applied a 
displacement of the atomic positions along each of the eight normal-mode unit vec- 
tors, £i t o with a prefactor ranging from -0.1 to 0.1 A. In Fig. [I] we plot as an example 
the resulting energies for the sixth h g mode. The pattern generated by V^ 9 (°) should 
be 1+2+2 (a state separated by two pairs of twofold-degenerate states). The small 
residual splittings of these twofold degeneracies, due to the cubic crystal field and 
higher-than-linear couplings, give an estimate of the accuracy of the method. By 
standard linear fitting and comparison with Eq. (|i~T|), we obtained directly the linear 
dimensionless coupling coefficients gh g i ah g i collected in Table 0. We determined the 
sign of ah g i by applying a distortion along £ i>1; and comparing the splitting of the 
HOMO with the eigenvalues of V /tfl (^(ia^). Following the same procedure we de- 
rived the couplings for the g g modes, by applying here g_i distortions. The resulting 
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Figure 1: The splitting of the HOMO degenerate level when the molecule is distorted 
according to the 6th mode h g (1250 cm -1 ). Points: DFT data. Lines: linear fits, 
giving the coupling parameters gh g % and a>h g 6 of Table |[ 

linear coupling coefficients g 9g i are also collected in Table ||. 

For convenience, we also report in Table || the amount of optimal JT distortion 
pertinent to each mode when the HOMO level is occupied by one electron/hole, 
and the corresponding energy lowering E s for both D 5 and D 3 minima. Note in 
particular the large coupling associated to the lowest h g mode, the corresponding 
distortion leading to an energy lowering practically equal to its quantum fkv. 

Adding up the JT energy gain of the individual modes, we estimate the total 
classical potential energy lowering in Cgo- The D 5 minima gain E s =71 meV, while 
the -D3 minima gain only E s =22 meV (the contribution of the a g modes - 2 meV - 
being included in both cases). It is therefore apparent that the Cg~ ion will choose, 
at least within linear coupling, the D 5 distortion. As was shown in (Manini and 
De Los Rios 2000), the possibility of a switch to a nondegenerate A u dynamical JT 
GS occurs, for large coupling strength g > 6, only under the condition that the D3 
minima are energetically lower or equal to the D§ minima. This settles finally the 
issue of the dynamical JT GS symmetry of this molecular ion: it is a regular Berry- 
phase vibronic state of symmetry H u , like the parent electronic state (Manini and De 
Los Rios 2000). No level crossing to a nondegenerate A u state is predicted to occur 
forC+. 

A static JT coupling resolves the degeneracy of the H u level: in a distorted config- 



9 



hu Ti 
_i 

cm 


coupl. 

eV/ao 


9ri 


OL T i 

deg 


q s (D 5 ) 

a 


a 


E S (D 5 ) 
meV 


E S (D 3 ) 
meV 


A 

meV 


a g 
500 


0.03 


0.059 


- 


0.00 


0.00 


0.0 


0.0 


0.1 


1511 


0.63 


0.274 




0.01 


0.01 


1 o 

1.8 


1 o 

1.8 


3.5 


9g 
483 


0.31 


0.757 


- 


- 


0.04 


0.0 


1.9 


8.6 


567 


0.05 


0.102 


- 


- 


0.00 


0.0 


0.0 


0.2 


772 


0.67 


0.800 


- 


- 


0.03 


0.0 


3.4 


15.3 


1111 


0.90 


0.624 


- 


- 


0.02 


0.0 


3.0 


13.4 


1322 


0.43 


0.228 


- 


- 


0.01 


0.0 


0.5 


2.1 


1519 


1.08 


0.467 


- 


- 


0.01 


0.0 


2.3 


10.3 


kg 

261 


0.50 


3.042 


-0.1 


0.27 


0.00 


30.0 


0.0 


75.0 


429 


0.43 


1.223 


30.1 


0.07 


0.03 


6.0 


1.1 


19.9 


718 


0.75 


0.995 


89.4 


0.00 


0.04 


0.0 


4.9 


22.0 


785 


0.67 


0.784 


-2.3 


0.04 


0.00 


6.0 


0.0 


15.0 


1119 


0.32 


0.221 


76.6 


0.00 


0.01 


0.0 


0.4 


1.7 


1275 


0.93 


0.519 


28.0 


0.02 


0.01 


3.3 


0.5 


10.7 


1456 


2.09 


0.962 


28.1 


0.03 


0.01 


13.0 


2.1 


41.7 


1588 


2.15 


0.869 


-31.1 


0.03 


0.01 


10.9 


2.2 


37.1 



Table 2: Computed mode eigenfrequencies and e-v linear coupling parameters of the 
H u HOMO in Ceo- The JT distortion magnitudes g s and the classical stabilization 
energies E s are tabulated for both D$ and D3 distortions and for one hole in the 
HOMO. The largest total JT energy gain is realized by the D§ distortions. We also 
show the contribution of each mode to the resistivity and superconductivity total 
coupling A = X/Ni(0) defined in Eq. fl2"U| ) [iVi(0) = density of states per spin per band 
at the Fermi level]. 

uration it is possible to distinguish individual levels within the HOMO, energy-wise. 
Figure |2| depicts the square modulus of the hole wavefunction for a D§ minimum. 
The probability density appears to be localized on an equatorial conjugated band, 
where the poles are the opposite pentagons centered around the -D5 axis we chose 
among the six possible ones. 

As a check, with the same method used above to calculate the hole-vibration cou- 
plings of the HOMO (Fig. [1]), we also computed the electron- vibration couplings of 
the LUMO and obtained the values in Table |3|. The total static JT potential energy 
lowering is 41 meV, of which 3 meV due to the a g modes, and 38 meV due to the 
h g modes. These values are generally in line with those calculated by previous au- 
thors (Antropov et al. 1993, Varma et al. 1991, Lannoo et al. 1991), although there 
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Figure 2: Electronic charge density distribution of the hole in Cg~ at the D 5 JT 
minimum along a go distortion of any h g mode. Out of the six equivalent ones, this 
particular minimum is defined by the D$ axis drawn as a bold "spit" piercing the ball 
through two pentagons. The contours are drawn at 30% of the maximum density. 

are some differences in the details. Error bars in theoretical determinations of JT 
couplings of fullerene have proven surprisingly large, possibly reflecting and ampli- 
fying errors in the vibrational eigenvectors. Of course as is well known, somewhat 
larger energy gains are obtained when the true dynamical JT problem, including a 
full quantum treatment of the vibrons is considered (Auerbach et al. 1994, Manini 
et al. 1994). We shall leave this calculation in Cgg" for future work. 

When more than one electron/hole occupies the HOMO of Ceo, larger distortions 
and more important energy gains are expected. In particular, the a g contribution for 
Cgo~ is simply: 

E«°(n)=n 2 E:°(l), (12) 

which can become as large as 179 meV for n = 10 holes. Consider now the nontrivial 
JT part of the coupling, that to g g and h g modes. For n = 2 electrons/holes (spin 
singlet configuration) in the HOMO orbital, the distortions simply become twice as 
large as in the n = 1 case, with JT energy gains which are four times larger. Even 
though n = 1, 2 electrons/holes take advantage only of the D§ stabilization energy, 
additional electrons/holes can benefit from the extra HOMO splitting induced by the 
h g {r = 2) plus g g coupling. To test this, we relaxed completely the molecular structure 
in the 64-dimensional space of the 8 h g and g g modes, and determined the minima of 
the total potential energy, filling the five BO sheets as drawn in the insets of Fig. |3|. 
For simplicity we only considered at this stage low-spin configurations, generally the 
most favored by JT. The resulting energy gains E^ a+ha {n) are reported in Fig. |3|. 
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Table 3: Computed mode eigenfrequencies and e-v linear coupling parameters for the 
T lu LUMO of C 60 . For each mode, we report The JT distortion magnitudes q s , the 
classical stabilization energies E s , and the contribution to A = A/iVi[0]. 

The energy lowering is maximum for n = 4 (and n = 6): it is as large as 401 meV, 
compared to a modest E^ 3 (n = 2, 4) = 153 meV in Cg~ . Similarly to what happens in 
the Ti u LUMO case (see Fig. |3| and Manini et al. 1994), the half-filled configuration 
n = 5 is slightly unfavorable (by JT energetics) with respect to the neighboring 
n = 4 and n = 6 states. We find that the contributions of the g g modes, strictly 
zero for n = 1,2, are small (~ 1 meV) but nonzero in the 3 < n < 7 configurations. 
This indicates that for such large fillings the many-modes JT system (Manini and 
Tosatti 1998) can take some advantage also of the "losing" g g + h g (r = 2) part of 
the coupling (favoring D 3 minima for n = 1,2). The insets of Fig. |3|, indicate that 
the largest displacement of a single level in the H u HOMO, about 0.28 eV, is realized 
for n = 2. Finally, we have particle-hole symmetry Ef 9+h9 (n) = E^ g+h9 (10 — n), at 
an opposite minimum distortion q s (n) = —^(10 — n). For the special case n = 5 this 
means that the configuration drawn in the inset and the one obtained reflecting the 
5 levels through zero give both the same optimal energy (at opposite distortions). 

Application of the above results to real CgQ 1 " must await the inclusion of electron- 
electron Coulomb repulsion. Coulomb interactions will generally compete with JT 
coupling and favor high-spin configurations, which, in turn, are generally less fa- 
vorable for JT. (For example for n = 2 the triplet configurations has a JT gain of 
100 meV instead of 277 meV for the singlet.) We shall return to a more detailed 
description of Cg^ ions in later work. 
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Figure 3: The static JT energy gain (excluding the a g modes contribution) of Cg^, 
as a function of the filling n of the H u molecular level. Electron-hole symmetrical 
values are obtained beyond half-filling n = 5. The corresponding energies for the T± u 
LUMO (Cgo~) are also plotted for comparison. Note the much larger gains (up to a 
factor 3) in the H u case. For each n, the insets show the corresponding split HOMO 
electronic configuration (all in the same scale of eV). 

4 Solid-State Transport 



In order to relate the previously calculated electron-vibron coupling constants to 
relevant solid state physical quantities, we start by re-deriving its contribution to 
the transport relaxation time. We are interested in particular to the features of an 
orbitally degenerate band as in the charged Ceo case. The Boltzmann equation in the 
presence of a uniform and static electric field E reads 



evkv ■ E 



dn°(e kfl ) _ f dn ku 
de kfl 



Of 



(13) 



coll 



being and the occupation number and energy at momentum k for orbital /x, 
and v kfl = de kfl /dk. 

Within the relaxation time approximation, dn k ll /dt = Sn k ^/T kfl so that the varia- 
tion with respect to equilibrium becomes 



, dn°{e kfl ) _ - 
on k)li = ^ eVk » ' 



(14) 
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and the conductivity is obtained in the form 

n 1 



a 



— E T kLiVkfM ■ V kfl 



3V 



dn°(e ktl ) 



(15) 



V being the volume, measured in cell units. 



As is usual in the treatment of different scattering mechanisms within Fermi's 
golden rule approximation, one has to sum the inverse of the corresponding relaxation 
times. The collision term due to the electron-vibron coupling is given, within Fermi's 
golden rule, by 



dn 



k,/x 



V / coll 

^k—q,u (1 ^k,fA, 

-riha (1 - n k 



2tt 1 
JiZV 



E E 



'g Ti huj 7 



q i,T,m,v \ 

N lU 5 - e *- q v - hoj Ti ) + (iv 



v$ m) M\ 2 { 



Tl 

q,m 



+ 1)5 (e kfl - e k ^ qv + hu T 



N 1 



l q , m d ( e k-qv ~ tkfi - h^ Ti ) + (N£ m + 1)5 (e fc _ ? „ - e feM + hu Ti ) } 



(16) 



Here iV™ m is the Bose-Einstein occupation number of phonons of crystal momentum 
q, symmetry r = a g ,g g ,h g , mode i, component m. n k ^ is the Fermi occupation 
of HOMO band component fi = —2, . . . , 2 at the given temperature and chemical 
potential. At equilibrium, the collision term is zero. We expand to first order in the 
deviation from equilibrium. 

First, let us consider the case T 3> MaXj^cUj) but still T <C Tp (the Fermi tem- 
perature). In this case 



q,m 



> 1, 



and 



dn 



dt 



coll 



27T 1 ^ ^ 

~W ZV ^ ^ 

a ^ V V i,T,m,v 



' g Ti fiuJ T 



2 k B T 

KUJn-4 



(Sn k ^ - 8n P;U ) [5 (e fcM - e pu - hu Ti ) + 5 (e kfl - e pv + ftuj Ti )]} . 



(17) 



Since T ^> Maxj(?ta;j), after inserting (|14]) into (|17D and using the Boltzmann 
equation, we find that 



Tku ft 2V V • 

K h l P i,T,m,u 

( _ v pv - fT 



' g Ti huj r 



2 k n T 



where n is a unit vector parallel to the electric field. 



[8 (e kfl - t pv - hu Ti ) + 5 {t kp - t pv + hu T i)} \ , (18) 
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If we approximately take = we find, since the electron-vibron matrices are 
symmetrical, 

\2" 



1 7T ^ / g Ti Tuo Ti \ Tr 



V 7 



ft. I 9. /i ?>/,)_, I 



1 / _ _ + s _ | 

where d is the orbital multiplicity, 5 for the HOMO-derived band of Ceo- Therefore, 
the relaxation rate in the linear T regime turns out to be roughly proportional to 

a 2 -huj • Tr f(V r ) 2 " 
A= V , (20) 

i,T,m 

related to the conventional dimensionless ep-phonon coupling A by 

A = iV x (e F ) A , (21) 

where N^ep) is the density of states per spin per band at the Fermi energy. 

By using Eq. ( p0|) and the calculated e-v coupling parameters of Table |], we find 
for holes in the H u HOMO-derived band 

A+ ~ 0.277 eV. (22) 

This value can be compared with that calculated similarly, using our couplings of 
Table |] for electrons in the T\ u band 

A_~ 0.197 eV. (23) 

We conclude therefore that holes have a stronger scattering with intra-molecular 
optical phonons than electrons, by about a factor 1.4. 

To confirm the accuracy of the values computed for the HOMO case, it is useful to 
compare those for the LUMO with previous similar calculations. It is then necessary 
to take into account the factor d = 3 which is usually incorporated in the total density 
of states (Lannoo et al. 1991). Our X/d = 66 meV compares well with previous similar 
calculations for the Ceo LUMO (Antropov et al. 1993, Varma et al. 1991, Lannoo et 
al. 1991), ranging from 49 to 68 meV. On the other hand, empirical estimates based 
on a fit to photoemission data tend to give a value as large as 147 meV (Gunnarsson 
et al. 1995). The origin of this discrepancy is presently unclear. 



5 Superconductivity 

The dimensionless electron-vibron coupling which governs the transport properties, A 



of Eqs. (|20| , |2lD , is not in point of principle coincident with the parameter A determin- 



ing the superconducting properties. The latter must be determined by solving the 
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Migdal-Eliashberg equation with the retarded interaction mediated by the vibrons 
plus the electron-electron Coulomb repulsion. However, a simple estimate of the or- 
der of magnitude of A can be obtained by taking the unretarded limit, and imposing 
a Debye cutoff to the electron energies. We assume the electronic band operators to 
be related to the molecular creation and annihilation operators through the unitary 
(orthogonal) transformation 

C-nku ^ nil (^0 djikm 

n being the band index, a the spin. Let us define a matrix W by 



9 T i nuJ ri ,rr(m)\ 2 ttat Sri^T 



- E fc (m) ) = E ^ E (24) 

as well as 



A y IJ.V 1 Z — t liv z. — , - 

r 5—1 



(25) 

The Bardeen-Cooper-Schrieffer (BCS) gap equation for 

&>k,n = ( c k,1,n C ~k,l,n) ) 



reads 



A fc ,„ = -^E >Vn fc , mp =^ tanh , (26) 

£v r, m " 'it try \ £ / 



p,m J -'p,m 

where 



E Pi m = V (e pm - Ho) + A p, m5 

/^o being the chemical potential. The critical temperature is obtained by solving the 
eigenvalue equation 



Snm^kp - E W nk ,mp-: T tanh ^ C ~ EIL - = 0. 

* V p,m \ € pm — H0\ \ Z ) 

In general, the BCS gap equation Q26D leads to interference between the Cooper pairs 
belonging to different bands. That, in turn, increases the critical temperature relative 
to a situation in which the pairs do not interfere. We may therefore foresee two 
opposite limits of strongly interfering and of non interfering pairs which, respectively, 
over and underestimate the effective coupling strength A. 

If we assume that averages over the Fermi surface do not depend on the band 
indices (interfering pairs, corresponding, for example, to the choice of Lannoo et al. 
1991), then we can replace W with W, and we find that the critical temperature is 
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determined by the maximum eigenvalue of the matrix W in Eq. 
assumption, the superconducting A is determined through 

A = iVi(e F )Max eigenva i ue (W) , 
For our H ® (a + g + h) e-ph coupling the matrices W r are: 



Under this 



(27) 
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(29) 



while W» is trivially the unit matrix. We note that, even though the matrices W r , 
for different r's, do not commute, the eigenvector (1, 1, 1, 1, l)/v^> is an eigenstate of 
each of them. Moreover, for each r this eigenstate provides the largest eigenvalue 



(30) 



(for any m and r). This means that the totally-symmetric paired state, delocalized 
over all the five H u orbitals or bands, is favored by the couplings to all modes. It also 
implies that, contrary to molecular JT, the couplings to all modes cooperate evenly to 
this superconducting state, and contribute additively to A. Totally equivalent (even 
if at first sight apparently different) results were derived for the d = 3 case (K 3 C 60 ) 
in Refs. (Lannoo et al. 1991, Rice et al. 1991). We note, however, that the claim 
that orbital degeneracy enhances the superconducting A through a factor d (Rice et 
al. 1991) is not really justified, as one must at the same time reduce the density of 
states from total to single-band, a factor 1/d smaller. We also note that 



s Tw T 

/ y ' ' mi 



y w t y v t v r 

i—imn ' ' mn L-imn y mn r nm 

d d 



E m [(V 
d 



T\2l 



Tr [(V 



T\2l 



(31) 
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which shows the identity of the A computed for superconductivity to the one obtained 
for transport in Eq. (f20|). 

If, in the opposite limit, the pairs did not interfere between different bands, the 
effective A would be reduced by a factor d. Although we cannot identify a physical 
situation corresponding to this limit, we can assume that a general case will be 
intermediate between the limits (interfering/not-interfering pairs). For simplicity, we 
will stick here to the interfering limit. 

In summary, is there any enhancement of superconductivity due to orbital degen- 
eracy? We can still identify one possible source for that, namely Coulomb pseudopo- 
tential. In fact, we note that, although a large A due to tunneling of the Cooper pairs 
between different orbitals/bands, can be seen as orbital degeneracy enhancing the 
effective A to the highest eigenvalue of W, there is no corresponding enhancement 
of the repulsive Coulomb pseudo-potential /x*, at least within the Migdal-Eliashberg 
theory. The reason is that the main contribution to the Coulomb pseudo-potential is 
a charge-charge repulsion which does not include tunneling processes between differ- 
ent bands, and being band-diagonal it does not get enhanced. In conclusion, in the 
above restricted sense, orbital degeneracy may in principle favor superconductivity. 

6 Discussion 

We have presented a density functional calculation of the linear coupling of holes/el- 
ectrons in the H u and T lu orbitals of the fullerene molecule to the intra-molecular 
vibrations. The coupling to holes is strongest for the h g modes, and among those 
to the lowest-frequency mode h g (l) around 270 cm -1 . The linear static Jahn Teller 
distortion predicted for Cq by these couplings corresponds to a D 5 distortion, with 
an energy gain of 71 meV. 

The corresponding dynamical JT state expected with the calculated coupling pa- 
rameters is a regular Berry-phase vibronic state of symmetry H u , like the parent 
electronic state (Manini and De Los Rios 2000). There is no level crossing to a non- 
degenerate A u state, as would have hypothetical^ been possible on pure symmetry 
grounds, had the D3 minima been the stable ones (Manini and De Los Rios 2000). 

In order to connect with important solid state properties including transport and 
superconductivity we have formulated a theory of the Boltzmann relaxation time, 
and of the BCS-type pairing, suitable for an orbitally degenerate multiband case 
with Jahn Teller coupling. This confirms that the same parameter A determines both 
transport and superconducting properties of the multiband degenerate solid. As we 
previously observed, not all the computed (Table 0) are small parameters, and thus 
weak-coupling BCS theory is strictly not applicable for Cg^. However, the overall A is 
still moderate. Our calculation neglects couplings to acoustic phonons and librations, 
which in principle should also contribute to e-ph scattering. In addition, similarly to 
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Lannoo et al.'s calculation (Lannoo et al. 1991), we assume that the dispersion of 
the HOMO band and of the optical phonons has a negligible effect on the integrated 
value of A. (This assumption was tested and proved correct in Ref. (Antropov et al. 
1993) for the couplings to the LUMO band.) With all these approximations, Eq. (|27|) 
should provide a semi-quantitative estimate of the total e-ph scattering. 

Assuming conservatively a total average density of states of ~ 10/0.6 eV ~ 17 stat- 
es eV _1 for the HOMO band (i.e. a single-band density of states iVi(O) = 1.7 stat- 
es eV -1 per band per spin), our calculated effective dimensionless A + for hole super- 
conductivity in C 60 is in conclusion about A + ~ 0.47. The Coulomb pseudopotential 
fi* is not available yet, but possibly in the same range of values as for negative Ceo 
[fi* ~ 0.2-=-0.3 (Gunnarsson et al. 1995, Gunnarsson and Zwicknagl 1992)]. With this 
value of A + , weak coupling would predict T c ~ 1.14 Tioj^ k^ 1 exp[— 1/(A+ — //*)] ~ 40 K 
for /i* = 0.2, and T c ~ 5 K for /i* = 0.3 (assuming a typical phonon energy u>d of 
about 1500 K). This seems of the correct order of magnitude, although somewhat 
on the low side, in comparison with T c = 52 K found experimentally. However, it is 
difficult to justify weak coupling in this case. 

The corresponding value A_ which we obtain for electrons in the T ltt orbitals is, 
assuming the same bandwidth of 0.6 eV, thus again iVx(0) = 1.7 states eV -1 per 
band per spin for the T lu band, A_ ~ 0.33. The factor A + /A_ = 1.4 of holes relative 
to electrons is in qualitative agreement with a larger T c of the former. Assuming 
the same typical phonon frequency, BCS would predict here T c ~ 0.8 K for fi* = 
0.2, and T c ~ K for fi* = 0.3. That is obviously way smaller than the observed 
superconducting T c = 10 K found experimentally in the field emission transistor (FET) 
experiment, let alone the higher values found in the fullerides. 

Coming to transport, the measured resistivities for holes are larger than those of 
electrons, and this also agrees with a larger A value. Calculation of the T-linear high 
temperature resistivity 

_ Af r T 
P ~ 4^5 

would however predict a moderately larger value for hole- than for electron-doped 
Ceo, at least assuming (somewhat arbitrarily) the same plasma frequencies for the 
same carrier densities. Quantitatively, Batlogg's FET data (Schon et al. 2000) differ 
strongly from this expectation. They, first of all, indicate a nonlinear temperature 
dependence, closer to T 2 ; secondly, they show values about 5 times larger for holes 
than for electrons. While there are second order processes (see Appendix) that would 
indeed yield a T 2 resistivity at low temperatures, we do not believe that they may 
explain the discrepancy here. Zettl and coworkers (Vareka and Zettl 1994) proposed 
that the apparent T 2 in the electron resistivity is an effect of thermal expansion, and 
showed that a linear T increase is recovered at constant volume, for negative doping. 

Recently Goldoni et al. measured by EELS the plasma frequency in K 3 C6o- They 
found it slowly decreasing with temperature, its width growing approximately quadrat- 
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ically with T. These data support the view that the T 2 resistivity is directly related 
with a T~ 2 decrease of relaxation time, most likely linked with lattice expansion. 
It seems plausible that a similar physics could apply to holes too, in which case 
the predicted constant-pressure relaxation-time drop with temperature would also be 
non-linear, and quantitatively larger than the electron case. If, on the other hand, it 
became possible to obtain the constant- volume inverse relaxation time and resistivity, 
then, assuming the same plasma frequency cu p , its increase should be linear with T 
with a slope 1.4 times larger than that of negatively- charged C 60 . This conjecture 
must await experimental test. 
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Appendix: Quadratic Resistivity 

Within Fermi's golden rule, the vibron contribution to the electrical resistivity is 
exponentially decreasing if temperature is much smaller than the vibron frequencies. 
However, the above result is not true any more if higher order corrections are taken 
into account. Indeed, at second order, the electron- vibron coupling generates an 
effective electron-electron interaction. Since the electrons involved lie on a shell of 
width T around the Fermi energy, the vibron-originated electron-electron interaction 

r,i,m 4 U ri ~ W 

^E E E^ m) («)^'(«)<^X/A<^^ , (32) 

kpq fivy/3 cr,a' 

acts as if it were effectively unretarded, u = 0. This interaction induces a T 2 inverse 
relaxation time at low temperatures. A rough estimation of the order of magnitude 
gives 

1 (k B T) 2 AT 3 _ f g 2 Tt huj Tl \ ( g 2 Tlj Tiu Tlj \ 
r K ~VT Nl{eF) .. \ A — 4 — J ^ 4 J 

z,j,t,t' ,m,m \ / \ / 

* Tr (v T ^V T ' < - m '^ 2 — Tr (v T ^V T '^ m '^V T ^V T '^ m '^ j 
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